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Abstract 

The results of an experimental exahiinaliOT of the 
vortex structures shed from a low profile “wishbone” gen- 
erator are presented. The vortex generator height rela- 
tive to the turbulent boundary layer was varied by testing 
two (Merently sized models. Measurements of the mean 
three-dimensional velocity field were conducted in cross- 
stream planes downstream of the vortex generators. In 
all cases a counter-rotating vortex pair was observed. In- 
dividual vortices were characterized by three descriptors 
derived from the velocity data; circulation, peak voriic- 
ity, and cross-stream location of peak voiticity. Measure- 
ments in the crossplane at two axial locations behind the 
smaller wishbone characterize the downstream develop- 
ment of the vortex pairs. A single region of streamwise 
velocity deficit is shared by both vortex cores. This is 
in contrast to conventional generators, where each core 
coincides with a regicm of velocity deficit 

The measured cross-stream velocities for each case 
are compared to an Oseen model with matching descrip- 
tors. The best comparison occurs with the data from the 
larger wishbone. 


Background and Research Objectives 


\brtex generators are used to provide a measure of 
flow control. Most often the goal is to achieve some 
resistance to boundary layer separation or to reduce the 
deleterious effects of locally strong secondary flows. 

In external flow situations, such as that encountered 
on airfoils, the most common application is upstream of 
flight control surfaces where boundary layer attachment 
is often critical to flight poformance. 

In internal flows vortex generators are used to pre- 
vent flow separation and reduce total pressure distortion. 
These effects occur often in inlet ducts and diffusers, 
due to such factors as duct centerline curvature and large 
changes in duct cross-sectional area. 
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The key to vortex generator performance is in the 
mixing and secondary velocity field created downstream 
by the shed vortex structures. If prc^rly situated in the 
flow field, the helical motion of the fluid in the vortex 
forces high energy fluid of the freestream into the slower 
moving fluid of the boundary layer while countering the 
tendency of naturally occurring secondary flows to distort 
the flow field by the opposite mechanism, namely to 
displace boundary layer fluid away from the wall and 
into core regions of the flow. 

Since the vortex generator presents an obstruction in 
the flow, there is a drag penalty paid for its use. Con- 
ventional blade type or delta wing type vortex gen^ton 
are found to work best when the vortex generator height 
above the flow surface, /i, is approximately that of the 
local boundary layer thickness, 6. Recent work in the de- 
sign and performance of vortex generators has produced 
low profile vortex generators that meet or exceed the mix- 
ing and strength perfcxmance of more conventional types. 
The height ratio, h/6y for a low profile vortex generator 
is typically between one and three-tenths, thereby pro- 
viding a substantial reduction in the drag penalty paid for 
its use. Low profile vortex generators have been used 
recently to improve the high angle of attack behavior 
of missiles,^ control a shock-boundary layer interaction,^ 
and improve the total pressure distortion and recovery 
performance of a diffusing S-ducL^ Figure 1 illustrates 
the use of a “wishbone” generator, which is a particu- 
larly effective low profile design. The secondary velocity 



Figure 1 The flow structures shed from a vortex 
generator are examined in a downstream crossplane* 

structure of the shed vortices is measured in a crossplane 
downstream of the vortex generator. As indicated in Fig- 
ure 1, this structure is a pair of counter-rotating vortices 
with the flow between vortices directed away from the 
wall (an “upflow” pair). This structure was verified in a 
flow visualization study conducted by Lin et al^ 
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Our qualitative description of embedded vortex struc- 
ture can be quantified using three parameters or “descrip- 
tors”. For any embedded vortex i these are: 

1 y.) — Location of the vortex core center. Its 

variation with axial position x represents the trajec- 
tory of the vcKtex. 

2 , — Peak vorticity. The peak vorticity is lo- 
cated at (zi, y,) and its magnitude is indicative of 
the vortex concentration. 

3. r, — Circulation of the v(Htex. Ti represents the 
“spinning strength” of the vortex. 

Each of these three descriptors are derivable from the sec- 
ondary velocity data taken in a downstream crossplane, 
as illustrated in Figure 1. A recent study conducted at 
NASA Lewis Research Center^ focused on the quantita- 
tive description of the vortices shed from conventional 
vortex generator blades. In this study we conduct a sim- 
ilar examination of the low profile wishbone generator. 
The objectives of this study are summarized as follows: 

1 . Mount a large low profile wishbone genaator in the 
thick boundary layer of NASA Lewis Research Cen- 
ter’s Icing Research Tunnel (IRT). At a suitable po- 
sition downstream of the vortex generator conduct 
highly resolved crossplane velocity surveys to char- 
acterize and quantify the vortex structures produced 
by the generator. Conduct these measurements on 
two different models, varying h/6. For one model, 
examine two downstream crossplanes to study the 
streamwise development of these structures. 

2. These measurements have the additional benefit of 
providing workers in computational fluid dynamics a 
basis for modeling low profile wishbone genaators 
in their studies. Simple mathematical models, devel- 
oped in terms of the vortex descriptors, are drived 
from the data to rqrresent the measured vortex struc- 
tures as velocity and vorticity distributions. These 
models are then compared against the data. 

Facilities and Procedures 
The Test Facility 

The investigation of tip vortices shed from vortex 
gen^tors generally involves a grid spacing having a 
scale similar to that of the vortex viscous core diameter. 
In the near- wake region of a low profile vortex generator 
this is roughly the same size as /i, and thus only a frac- 
tion of the local boundary layw* thickness. Typical test 
section turbulent boundary layers gen^ted in facilities 
at NAS A-Lewis are on the order of a centimetCT or so in 
thickness. To obtain the fine scale grid resolution desired 
for this study it was necessary to locate a test facility ca- 
pable of generating a nominally two-dimensional turbu- 
lent boundary layer with 6 in excess of 7 or 8 centimeters. 
This requirement was determined as a result of balancing 



Figure 2 A schematic illustration of 
NASA Lewis’ Icing Research IhnneL 

such factors as time available for testing and the size of 
five-hole probe tips available for use within Lewis’ Inlets, 
Ducts, and Nozzles workgroup. NASA Lewis Research 
Center’s IRT is a suitable candidate, having a turbulent 
test section wall boundary layer of the required thick- 
ness. The IRT is a large subsonic test facility designed 
and equipped to support the testing of low speed mod- 
els. It is a closed-return atmospheric-type tunnel with a 
rectangular cross-section at every axial station. Figure 2 
is a plan view of the IRT showing its various operational 
elements. In this study tunnel total pressure and to- 
tal temperature are ^proximately atmospheric. The test 
section dimensions are 1.8 by 2.7 meters with velocities 
ranging between 20 and 140 metcrs/sec. Corresponding 
mass flow rates are between 100 and 500 kilograms/sec. 

Operation of the wind tunnel is fully automatic. 
Pressure and temperature instrumentation located in the 
test section provides the data input to aid in setting and 
controlling tunnel flow conditions. Each flow measuring 
device, whether operational or research, is monitored by 
a central minicomputer based operating system. 

More information concenung the design, op^ations, 
and capabilities of the IRT may be found in the rqx)rt of 
R. So^er and C. Andracchio.^ 

Test Parameters and Research Instrumentation 

Figure 3 is an illustration of the test section showing 
the coordinate system in use. The cartesian coordinate 
axes are centered on the tunnel tc^ wall at an axial loca- 
tion corresponding to the boundary between the contrac- 
tion and test section. The model vortex genmtor, probe, 
and probe traversing hardware are mounted on a hatch in 
the test section ceiling. The leading edge of this hatch 
is located at x » 3 meters. This test was conducted 
at a freestream Mach number of 0.2. The test section 
freestream velocity, Uoo* is 2 q>proximately 70 meters/sec. 

Boundary layCT surveys conducted in the empty test 
section indicate a nominally two-dimensional turbulent 
profile over the central two-thirds span of the top wall. At 
the axial location corresponding to the mounting position 
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Figure 3 The IRT test section showing the coordinate 
system used and the location of the test hardware. 


of the vortex gen»ator model <5 « 9 centimeters, and 
the Reynolds numbw is Re = Uoox/v « 13 x 10®. 
No significant spanwise variations in profile shape or 
thickness woe detected in the center span region of the 
test section ceiling. 

The vortex genmttor models tested in this study are 
patterned after the Wheeler wishbone design.’ Figure 4 
is an illustration of the geometry. Two models were 
constructed from a block of aluminum alloy using a wire 
cutting (EDM) process. The dimensions of both models 
are given in Figure 4. 

Figure 5 illustrates the arrangement of hatch, vortex 
graerator model, velocity probe, and traversing bearing. 
A two axis (z, y) mechanism located above the hatch 
moves the probe stem in the crossplane survey grid, 
sealed by the firee-floating “double-bearing” rig. 

The velocity probe is a five-hole Pitot probe with 
a tip diametCT of 0.23 cms. This probe tip is mounted 
off-axis on a goose-necked stem to a position upstream 
of the circular bearing seal line. The probe is used in 
a non-nulling capacity and is calibrated to determine all 
three components of the mean flow velocity following the 
procedure outlined by Treaster and Yokum.* 



Figure 4 Small wishbone generator geometry. 
The length dimensions of the larger generator 
are 2.5 times the values listed here. 



Figure 5 The arrangement of model generator 
and instrumentation on the test section ceiling. 

Right and left side flow symmetry in the crossplane 
is expected in the wake of the vortex genoator. To 
taVe. advantage of this, the survey grids examined in this 
study are not centered with respect to the model, but 
are positioned to favor the right half. The large vor- 
tex gen^tor model is mounted with its trailing edge at 
{x, z) = (3.24, 0.00) m. The 15.1 x 13.5 cm (z, y) 
crossplane survey grid used for the larger model is cen- 
tred at (x, z) = (3.47, 0.04) m. This is one chord length 
downstream of the model. 

The smaller vortex genmtor is mounted at two 
I locations. The (x, z) trailing edge positions are at 
(3.10, 0.00) m and (3.38, 0.00) m. The axial variation 
of S between these two x positions is negligible. The 
1 1 .4 X 10.8 cm grid used for the smaUw model is centered 
at (x, z) = (3.47, 0.04) m for survey locations of four 
chord lengths and (me chord length downstream of the 
mcxlel, respectively. The grid spacing for both models is 
Az = Ay = 0.64 cms. 


Experimental Results 


Small Vortex Generator 

The crossplane vel(x:ity field measured one chord 
length downstream of the small vortex generator 
(h/8 — 0.30) is shown in Figure 6a. A concentrated pair 
of upflow vortices is evident in the vector plot. 

The three structural descriptors of the RHS vortex 
are derived firom the measured crossplane velccity field il- 
lustrated in Figure 6a. The vel(x:ity field is first converted 
to a streamwise v(Mticity field following the relation: 



Finite difference formulas are used to rqmesent the spatial 
derivatives in Eq. (1). Figure fib is the vorticity plot 
correspontfing to the veltxity field depicted in Figure 6a. 
^jmax located at a grid point having coordinates (z, , yj). 
For the RHS vortex illustrated in Figures 6a-b: 

0)7*“ = 7586 sec- \ (z<, y<) = 1.44, 2.56 cms. (2) 
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The center of the vortex appears to be located on the right 
shoulder of the vortex genmitor. The circulation, Pj , is 
calculated by first isolating the region of core vorticity in 
the field. This is done by referring to the contour 
plot in Figure 6b. A path enclosing the region of core 
vorticity is defined. The circulation is then calculated 
according to: 

Ti=^Vds, (3) 

path 


whoe is the velocity vector in the crossplane, and 
s refers to the path coordinate. By using closed paths 



Figure 6a - Transverse velocity. 



Figure 6b - Stream wise vorticity. 


y cBis 



Figure 6c - Primary velocity. 


composed of line segments in the z ot y coordinate di- 
rections, the circulation is easily detomined. The circu- 
lation of the vortex illustrated in Figure 6b is found to 
be r< = 1.561 m^/sec. Contours of streamwise veloc- 
ity ratio are plotted in Figure 6c. The boundary layer 
is slightly thickened in the vicinity of the embedded up- 
flow pair, as would be expected. An unexpected result is 
the lack of a streamwise velocity deficit coinciding with 
the location of each vortex center. A single large region 
of streamwise velocity deficit encloses both vortex cores 
and is centered at a crossplane location coincident with 
the tip location of the vortex gen^tor. 

The ccMTespotuling results recorded four chOTd 
lengths downstream of the small vortex generator are 
shown in Figures 7a-c. The descriptors of the RHS vor- 
tex at this location are: 

^mor _ sec"^, (zi, jti) = 2.08,5.12 cms, 
r.- = 1.047 mVsec. 

Figures 7a and 7c illustrate the displacement of the 
vortex cores away from the wall, convected there by the 
strong secondary flows occurring in the upwash region 
of the vortex pair. As discussed in Wendt ei al.? the 
33% reduction in circulation may be attributed, in part, 
to the crossplane component of wall shear stress which 
applies a torque opposing the rotatirm of the vortex. An 
additional mechanism fw the obsoved circulation decay 
is the diffusion of vorticity between the closely spaced 
vortex cores. A steep gradient in streamwise vorticity 
between the LHS and RHS cores is evident in Figures 
6b and 7b. In addition to reducing circulation and peak 
vorticity, diffusion of the core vcxticity has significantly 
increased the core size as can be seen by comparing 
Figures 6 and 7. Located under the core in Figure 7b is a 
r^on of “secondary” vorticity. This vorticity, oiqx>site 
in sign to that of the core or “primary” vcaticity, is 
induced by the viscous interaction between the secondary 
flows of the vortex and the walL Figure 7c again shows a 
single region of streamwise velocity deficit This region 
is now centered between, and above, the vortex cores. 

Large Vortex Generator 

The results obtained one clKH'd length downstream of 
the large vortex genoator {k/6 = 0.70) are illustrated in 
Figures 8a-c. The descriptors of the RHS vortex are: 

^max _ sec“* , (zj, yi) = 2.74, 8.15 cms, 
r, = 5.173 m^/sec. 


Modeling 

The large numbo’ of parameters to consider when 
designing a vortex generatcH' array for an aircraft com- 
ponent such as a wing or inlet has the implication that 




Figure 6 The results one chord length 
downstream of the small vortex generator. 
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experimental work on optimum designs is often slow and 
expensive. This fact has motivated a few workers in 
computational fluid mechanics to assist in the optimizing 
problem by including a means of representing vortex gen- 
erators in their codes. A simple and effective means of 
doing this is to employ a model for the crossplane veloc- 
ity or vorticity induced by the gen^tors. This is the 3p- 
proach taken in recent work by Anderson et who 

examined multiple vortex generator array geometries in a 
diffusing S-duct inlet using a parabolized Navier-Stokes 
(RNS) solver. A similar inclusicm of embedded vortices 
in a full Navier-Stokes (FNS) code was implemented by 
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Figure 7 The results obtained four chord lengths 
downstream of the small vortex generator* 
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Figure 7a - Transverse velocity. 



Figure 7b • Streamwise vorticity. 



Figure 7c • Primary velocity. 
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Figure 8b - Streamwise vorticity. 
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Figure 8c - Primary velocity. 

Figure 8 The results obtained one chord length 
downstream of the large vortex generator. 

Cho and Greber^^ for a constant area circular duct and 
a diffusing S-duct geometry. The advantage of this ap- 
proach lies in the fact that a newly formed vortex may 
be accurately represented in this manner. 

Experimental workers concerned with embedded 
vortices have often noticed the close similarity of 
observed crossplane vortex structure to simple two- 
dimensional models of vortices. Inviscid or “potential” 
models were first applied to results obtained in subsonic 
diffusers by Taylor and Grose. Later, Pearcy^"^ would 
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develop embedded vortex interaction models based on 
the inviscid rqjresenlation of embedded vortex structure. 
Eibeck and Eattm’* compared the structure of a single em- 
bedded vortex to that of the patchwork “Rankine” model 
vortex. Studies by Pauley and Eaton** and Wendt et alP 
have made comparisons to the “ideal viscous” or “Os- 
een” model. 

Lets now examine the equations of the Oseen model 
and show the comparison between the model and the 
present results. 

The Oseen Model of the Embedded Vortex 

The two-dimoisional Oseen model of a viscous vor- 
tex rqHesents the time dq)endent decay of a potential 
vortex where the velocity at the origin (r = 0) is forced 
to zero at time t = 0. A single isolated vortex cen- 
tCTed on the crossplane origin has velocity components 
(in cylindrical coordinates); 


constructed. Eqs. (8) represent the secondary velocities 



Figure 9 The constructioD of an Oseen 
model of a single embedded vortex. 


vt 


Vr = 0 , 





( 6 ) 


where r is the vortex circulation and u the laminar 
coefficient of kinematic viscosity. Following Squire** 
the unsteady solution is transfexmed to a steady one by 
relating the decay time to the distance between the vortex 
generator tip and the crossplane of interest {x - xo): 


X - Xq 


(7) 


For an isolated turbulent vortex i located at (z,, yi) the 
velocity components in the crossplane can be written in 
terms of the vtHtex descriptors.*’ In cartesian coordi- 


nates: 


Vi = 
Wi = 


Tijz-zi) ^ 

r<(y-yQ p 


( 8 ) 


where: 

Ri ={z- Zif + {y-yif , (9) 


and: 

.max 

The vorticity field is given by: 


Vi and for the vortex above the wall. The equations for 
the image vortex are identical excq)t for the sign switch 
on r, and y, . Denote these velocities v,mi and iy,mi • 
The equations representing the embedded vortex are then 
simply: 

= + 

W = Wi+ Wim 1 • 

Figure 10 illustrates the construction of an embedded 
array of N vortices. The secondary velocities of this 



Figure 10 The construction of an Oseen 
model of an array of N embedded vortkies. 

flow field are: 


wf = 


(1 - Fi). 


( 11 ) 


V = t?i + Vim 1 + ^2 + Vim 2 + • * ' 
h + Vim JV, 

K; = tVi -f Wim 1 -I- IV2 + 2 “I 

1- Wpf -j- Wim N • 


(13) 


We can superimpose solutions for a representation 
of the wall and neighboring vortices. Figure 9 illustrates 
how a representation of a single embedded vortex is 


6 





Model and Data Comparisons 


y cms 


The model is constructed by matching the vortex 
descriptors in the equations to those previously listed for 
the data (Eqs. 2, 4, and 5). The circulation and peak 
voiticity descriptors for the LHS vortex in each case are 
chosen to be equal (but opposite in sign) to the descriptor 
values determined for the RHS vortex, (z, , y,) for the 
LHS vortex is found by symmetry with respect to the 
line (x, y = 0, z = 0). 

Figure 11a illustrates the crossplane velocity field 
of the model when the vortex descriptors are matched 
to those obtained from the data taken one chord length 
downstream of the small vortex genCTator, This can be 



Figure 11 The model velocity field (11a) one 
chord length downstream of the small vortex 
generator, and a contour plot of the variance 
between the measured velocities and the model (11b). 



Figure 12 The model velocity field (12a) four 
chord lengths downstream of the small vortex 
generator, and a contour plot of the variance 
between the measured velocities and the model (12b). 

small vortex generator. Figure 12a may be compared to 
Figure 7a. The best comparison here occurs just below 
the vortex core, and to either side. 

Figures 13a-b iUustrate the model and variance ob- 
tained one chord length downstream of the large vortex 
gen^tor. Figure 13a may be compared to Figure 8a. 
The model provides a better match to the data in this 
case, as evidenced by the large regions of the flow field 
(in the vortex upwash, downwash, and above the core) 
where the variance is less than 15%. 


Summary 


compared to Figure 6a. One point to note is that the 
near wall velocities arc a good deal weaker in the model 
rqpresentation of the data. Figure lib is a contour plot of 
the variance between the data and the model defined as: 


variance = 


V<iata ^model 

I Vdatall 


(14) 


The best comparison occurs just above the RHS core, 
where the variance is less than 15%. 

Figures 12a-b illustrate the corresponding results for 
the crossplane field four chord lengths downstream of the 


The low profile wishbone generator sheds a strong 
pair of counter-rotating vortices. Measurements of ve- 
locity in a downstream crossplane form the basis for 
characterization of the vortex structures. Vortex struc- 
ture was quantified by three descriptors. Vortex strength 
was charaetCTized by its circulation or integrated stream- 
wise vorticity, vortex concentration by the magnitude of 
peak vorticity, and the vortex trajectory by the location 
of peak vorticity. 

Measurements made at two axial locations illustrate 
the streamwise development of the shed vortices. The 
counter-rotating pair of vortices was obs^ed to lift off 
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Figure 13 The model velocity field (13a) one 
chord length downstream of the large vortex 
generator, and a contour plot of the variance 
between the measured velocities and the model (13b). 

the surface quickly. A single region of streamwise veloc- 
ity deficit was located between vortex cores. Although 
the observed three-dimensional flow structure was com- 
plex, the flow in the crossplane is well rqiresented by the 
two-dimensi(xial Oseen model, particularly in the case of 
the larger vtxtex generator. 
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